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It is shown that in magnetic superconductors with spiral magnetic order 
the Josephson current has an additional contribution which depends: (i) on 
the relative orientation (magnetic phase) 9 = 9l — 9r oi magnetizations on 
the left (L) and right (R) banks of the contact, (ii) on the junction helicity 
X = XlXr, (with spiral helicity XLiR) = ±1) , i-e- J = [Jc - cos 6^] sin 93 
with ip = ipi — ipf{. The ratio = J^/ Jc is calculated as a function of the 
superconducting order parameter A, the exchange field energy h and the 
wave vector Q of the spiral magnetic structure. The vr-Josephson contact 
can be realized in such a system in some region of parameters. Some possible 
consequences of this new phase relation is also analyzed. 
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Introduction - The physics of magnetic superconductors (MSG) is interesting due 
to competition of magnetic order and singlet superconductivity (SC) in bulk materials. 
The problem of their coexistence was first studied theoretically by V. L. Ginzburg 
in 1956, while the experimental progress begun after the discovery of ternary rare earth 
(RE) compounds (RE)Rh4B4 and (RE)Mo6X8 (X=S,Se) [0] with regular distribution of 
localized RE magnetic ions. It turned out that in many of these systems SC, with the 
critical temperature T^, coexists rather easily with antiferromagnetic (AF) order, with the 
critical temperature Taf, where usually one has Taf < 0- Due to their antagonistic 
characters singlet SC and ferromagnetic (FM) order can not coexist in bulk samples with 
realistic physical parameters, but under certain conditions the FM order, in the presence 
of SC, is transformed into a spiral or domain-like structure - depending on the anisotropy 
of the system 0], P]. These two orderings coexist in a limited temperature interval 
Tc2 < T < Tm (reentrant behavior in ErRh4B4 and HoMoeSs) or even down to T = Q K 
(in HoMoeSeg). Note, in ErRh4B4 one has = 8.7 K, T„ ^ 0.8 K, T^s ^ 0.7 K, while 
for HoMoeSg one has = 1.8 K, Tm ~ 0.74 K, T^a ^ 0.7 K - see Refs. §, §, §. In 
most new quaternary compounds (RE)Ni2B2C the AF order and SC coexist |^, while in 
HoNi2B2C, an additional oscillatory magnetic order exists competing strongly with SC 
leading to reentrant behavior p. Recently, the coexistence of SC and nuclear magnetic 
order was observed in Auln2 [0, where Tc = 0,207 K and Tm = 35 fiK, which was 
explained in terms of spiral or domain-like magnetic structure [Q. There are evidences 
for the coexistence of the AF (or FM) order and SC {Taf = 137 K, T^ < 45 K) in layered 
perovskite superconductor RuSr2GdCu208 0, where AF and SC orders are spatially 
separated in Ru-0 and Cu-0 planes, respectively. 

In all of the above MSC systems the exchange(EX) interaction between localized 
magnetic moments (LM's) and SC is much larger (the EX model) than the electromagnetic 
(EM) one. The latter is due to the orbital effect of the magnetic induction B ^ 47rM on 
SC electrons. 

In the following we study the Josephson effect in MSC with spiral magnetic order in 
the framework of the EX model. It is shown below that the Josephson current depends 
additionally on the magnetic phase 6 = 6l — Or of the magnetic order parameters, as 
well as on the spiral helicities Xl{R) = ±1 on the left and right surface. 

The microscopic theory of MSC [|| takes into account the interaction between LM's 
and conduction electrons which goes via: a) the direct EX interaction and b) the EM 
interaction where the dipolar magnetic field Bm(r) = 47rM(r) created by LM's acts on 
the orbital motion of electrons. The Hamiltonian of the system has the form 

H = Ha + Hbcs + ( iPr[i,\v)V„(i)i,(v)+ 



+ + (1) 

i 

Hq = Hq{p — |A) describes the motion of quasiparticles in the magnetic field B(r) = 
curlA(r), which is due to LM's and screening superconducting current. Hbcs = 
^/'bcs(^(i')) is the BCS pairing Hamiltonian with the SC order parameter A(r), Vex{r) = 
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h(r)cr is the EX potential where h(r) = J^i'^exi^' — ^i){g ~ l)Jj- Here, Jex(r) is the 
exchange integral between electronic spins a and LM's and J, is the total angular mo- 
mentum operator of the LM at the i-th site and a = (ai, a2, cts) are Pauli spin matrices. 
In case of magnetic anisotropy the crystal-field term HcF^^i) should be added to H. 

In the following a clean s—wave MSG is considered with spiral magnetic order with the 
wave vector Q along the 2;-axis, Q = QzZ = ±Qz, and the spiral helicity x = Qz/Q = ±1- 
The LM's are assumed to lie in the xy-plane due to the easy-plane magnetic anisotropy 
and the mean-field EX potential V^(r) =< Vex{r) > reads 

^y^)- yhe'^xQz+e) or 

Here, h = nm{g — l)JexiO) <| J |> and Jex{0) is the q = Fourier component of Jex(r — rj), 
rim is the concentration of regularly distributed LM's - see ||^. For further purposes we 
define ho = he^^ where the magnetic phase 9 characterizes the orientation of the EX field 
at the surface z = 0. The Josephson current between two MSG, both with spiral magnetic 
order, depends on anomalous Green's functions Fl^^,^{xi,X2) =< Ttpl_^{xi)ip fo-a (^2) > 
with X = (r, r) and a =t, |, which are calculated in [i] 



Fl.ik, k'; = 5(k - kOA "" + ^' (3) 

'^n + ^IkJ K + ^2kJ 



Fl,ik + Q, k'; c.) = -6{k - k02A/.g f| (4) 



The excitation spectrum -Ei,2k is given by Ef2k = ^k + Pk + I ^ 

T2V4(Pk + ^')+ I A 1'^', where = [^(k + Q/2) + ^(k - Q/2)]/2 and pk = [^(k 



Q/2) — ^(k — Q/2)]/2, i.e. ^k = ^(k) and Pk = ■ Q- The pair function Fj, is obtained 



from fJ^ by replacing —>■ hg, i.e. 6 —>■ —6 and Q — Q . In the following we assume 
that Q <^ kp. Note that in the MSG system with spiral magnetic order triplet pairs with 
amplitudes Fj^ and (which depend on he) live short time inspite of the absence of 
triplet pairing. Below is demonstrated that the finite life-time (of the order h/ A) of the 
triplet pairs leads to nontrivial effects in the Josephson current. 

The above theory predicts that if the EX interaction dominates then SG modify 

the FM order into the spiral with the wave vector Q, ^ Q kp, and with zero average 
EX field < h >^o~ (^0 is the coherence length). The two ordering coexist in a narrow 
temperature interval in ErRh4B4, HoMogSs, or up to T = i^' in HoMoeSes. In Auln2 
SG coexists with the nuclear (spiral) magnetic order up to T = -ft^ 0, [§. 

Josephson current - We consider tunneling contact (MSG-I-MSG) between two (left- 
L and right-i?) MSG both with spiral magnetic ordering with the wave vector Q,l,r, 
the SG order parameter Al,r =\ Al,r \ e*'^-^-« and the EX field he^^^^ = /iL(ij)e*^^(«\ 
respectively - see Fig.l. For simplicity it is assumed that: (i) | Al \ = \ A/j | = | A |, 
hi = hn = h while Lp = Lpi - '^r ^ ^ and 6 = 6l - Or 0; (ii) | Ql | = | Qr \= Q 
where Ql,h = Xl,rQz are orthogonal to the tunneling barrier, where helicities of MSG 
are Xl(r) = ±1 for Ql,_r along (+) or opposite(-) to the z-axis. In that case one has 
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Pk,L(R) = ■ Ql{r) = Xl(r)P = Xl(R)Po COS f3, where (3 is the angle between vp and 
2;-axis. 
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FIGURE 1. The Josephson contact / between two magnetic superconductors Sl and Sr with spiral 
magnetic orders. The corresponding exchange fields /il.a at the surface make angles 9l,r with the j/-axis. 
The wave vectors Ql,r of the spirals are along the z-axis 



The tunnehng of a left-side electron with momentum and spin (ki, a) into a right-side 
one (kij, 0") is described by the standard tunneling Hamiltonian ||10|. For simplicity it is 
assumed, as usual, that the tunneling amplitude Tk^^k^ is weakly energy and momentum 
dependent, i.e. Tk^,kj{ ~ T^QiJiL^zkR^z)- The Heaviside function takes into account that 
before tunneling the left electron moves toward the barrier, while after it moves as the 
right electron away from the barrier, and vice versa. The assumed Tk^^k^ is more suitable 
for diffusive barrier (incoherent tunneling). The extension of the theory to the momentum 
dependence Tk^^k^ is straightforward. 



The standard theory of the Josephson effect []10l applied to MSG systems gives the 



Josephson current J(v9, 9) and the energy of the junction Ej{ip, 9) 

J{ip, 9) = [Jc — Jx cos 9] sin ip, (5) 

Ej{(p, 9) = [1 — cos 9] cos Lp + const, (6) 

where, ip = ipi — V^_r, 9 = 9l — 9r, $0 is the flux quantum and Ry. = Jc/ Jx- The 
first term in the bracket (~ Jc) is the standard one due to the singlet pair tunneling 
Zlki,kfl,t^„ I TkL^R P -^u,L(kL,t^n)-^Ti,fl(kR, -uJn)- In the calculation the summation over 
k is replaced by the integration over ^ and p(= po cos/5 = poy), i.e. 

Ar(o) rpo 



2po 



dp / da-) (7) 
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iV(0) is the density of states on the Fermi leveL After the integration over ^ one obtains 
Jc (in the following | A |= A) 

oo „i 

J, = AeTi^N^Q) I To p A^T V[ / I{u;^,y)dyf 

n=l 



I{ujn,y) = 



at 



ttn = uj^ + ply"^ + h'^ + A^, po = Qvp, is the Fermi velocity and T is the temperature. 

The second term in -Eg.® (~ J^) depends on the relative magnetic phase 6 = 6l — 6r 
of the EX fields at the barrier surfaces, and on the Junction helicity x(= XlXr) = 
±1. Both are due to the tunneling of Cooper pairs being short time in the triplet 
state. is due to the term Eki,kfl,c^„ I ^k^.k^ P [T|^_^(ki, a;„)[T|^_^(kjj, -u;^)]* + 

^(ki, ci;„)[-F^ ^(kij, — ci;„)]*]. After the ^-integration reads 

oo „i 

= lQen''N^{0)tlA^h^Tj2{[^n / K{iUr.,y)dyf+ 

n=0 



+XlXb[Po / yK{ujn,y)dyf}, 
Jo 



i^K, y) = . (9) 



From follows that = for h = 0, while the standard Josephson current 

Jc is finite at /i = reaching its maximum. > for the total Junction helicity 
x{= XlXr) = 1) while for x = — 1 it can be negative depending on A, h, po - see discussion 
and Figs. 2b, 3b below. In order to calculate Jc and equilibrium values of A, h, Q, 
which minimize the free-energy F{A,h,Q), are needed. As it is shown in 0] F{A,h,Q) 
(and A, h, Q) depends on microscopic parameters kp,VF,nm, Jex, ^o, lattice structure, 
etc. what shall be not studied here. 

In the following the ratio R^{m,p) = Jx/Jc is calculated numerically as a function of 
m = A/h andp = (po//i)(= Qvp/h) at the temperature T = 0.1 A. m and p are supposed 
to be equilibrium values. The results are shown in Fig.2 for p = const 
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FIGURE 2. The ratio R^{m = const, p) = J^/ Jc for various m = 1; 0.1; 0.02. (a) x = +1; (b) X = -1- 
For I |> 1 the tt - contact is reahzed. 



and in Fig. 3 for m = const. 




FIGURE 3. The ratio R^{m,p = const) = J^/ Jc for various p = 2; 1; 0.2. (a) x = +1; (b) X = -1- For 
I i?^ |> 1 the TT - contact is reahzed. 

Discussion 

(i) n-junction: - The above theory predicts that, besides | i?^ |< 1, the case | i?;^ |> 1 
can be reahzed depending on m and p - see Figs. {2 — 3). From comes out 

that for I |> 1 the vr-junction is reahzed, i.e. min i?j is reached for = vr. ^ 
imply also that for x = — 1 the case < — 1 is realized in some region of m and p 
(see Figs. (2-3)) and the vr-contact is realized for those 9 with cos9 < — 1/ | R^ \. For 
Ry^ > 1 the TT-contact is realized for cos^ > V I -^x I- '^^^ favorable range of parameters 
m and p for the realization of vr-junction are seen from Figs. [2 — 3) and it lies in the 
region m < 1, p < 2. Note that m(= A/hex) is a measure of the strength of the SC 
order parameter with respect to the magnetic order parameter (exchange energy), while 
p(= Qvp/hex) = lex/Lspirai measurcs the inverse of spiral period {Lgpirai = 2ti/Q). This 
means that MSG systems with spiral magnetic order which are more ferromagnetic-like 
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are favorable for the realization of the vr-contact. On the other side the latter property 
is less favorable for the coexistence of SC and spiral magnetic order. In that respect it is 
worth of mentioning that in ferromagnetic SC ErRh4B4, HoMogSs, HoMoeSes and Auln2, 
m varies from m > 1 to m <C 1 by lowering T from Tm, while p ^ 10 is realized thus 
making | \ small in these systems 

The proposed theory holds also for MSG with AF magnetic order (the limiting case 
of the spiral order) with an easy-plane magnetic anisotropy. For instance in some AF 
heavy-fermion superconductors, like URu2Si2 with Taf ^ 17 K and Tc ~ 1.5 K, one has 
p ?a 1 — 2, m ~ (0.01 — 0.03), where small value of p is due to the small Fermi velocity [|lT 



If one assumes that in this system (anisotropic) s-wave SC is realized in absence of the AF 
order, then the latter changes SC in the way described above giving rise to short living 



triplet pairs, gapless superconductivity, power low behavior, etc. [|12[- Since in URu2Si2 
one has p ^ 1 — 2, m ^ (0.01 — 0.03) and | |> 1 this means that it is favorable for 
making vr-j unctions. In that respect other heavy fermions, like UPd2Al3 {Taf ~ 14 K, 
Tc^2 K) and UNi2Al3 {Taf ~ 4.6 K,Tc^ 1 K), might belong to the class described by 
the above theory [jl2[. UPta {Taf 5 K, T^ ^ 0.5 K) is also a candidate for such a vr- 



junction, if it can be described by the above theory, because p ^ 1 — 2, m ^ (0.01 — 0.03) 
in it. However, various experiments in UPta are well described by the unconventional 
superconductivity [|r^ - the type of order parameter is still under debate. 

(ii) magnetic phase 6: In the case of an easy-plane (x-y plane) magnetic anisotropy the 
magnetic phase difference 6{= 6^ — Or) can be tuned in the range (vr, — vr), for instance, 
by applying magnetic field on lateral surfaces of the left and right superconducting banks. 
So if the parameters {m,p) of the system allows that | R^ \> 1 then by rotating the 
external magnetic field one can tune the Josephson junction from the 0- to vr-junction. 
Note that contrary to the vr-junction based on d-wave high-Tc-superconductors junctions 
based on MSC system can be varied from 0- to 7r-junction by simple changing orientation 
of external magnetic field. Even in the case | R^ \< 1, when only 0-junction can be 
realized, by tuning 6 one can make significant changes of the Josephson current. 

The proposed MSC-I-MSC junction opens new physical possibilities if spatial and time 
variations of 6{x, y, t) are realized. It is to expect that magnetic and electric fields, applied 
in the contact, will change both phases 9{x,y,t) and ip{x,y,t). The elaboration of these 
effects in a form of coupled equations for ip{x,y,t) (a modified Ferrell-Prange equation) 
and for 6{x,y,t) is the matter of future researches. 

In conclusion, we demonstrate that in magnetic superconductors with spiral magnetic 
order the Josephson current depends on new degrees of freedom: (1) the relative magnetic 
phase (orientation) 6 = 6l — Or of the exchange fields (magnetizations) on the barrier, 
and (2) on the helicities of the left and right spiral, i.e. on x{= XlXr) = il- some 
range of parameters A, hex and Qvr the vr-junction can be realized by tuning 6 in external 
magnetic field. Even in the case of the 0-junction the Josephson current can be varied 
significantly by tuning 6 thus giving rise to interesting physics. Time and spatial changes 
of the magnetic phase 6 and the Josephson phase can be of potential interest for small- 
scale applications, of course if the MSC-I-MSC junction is realizable. 
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